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The geometry of electronic bands in a solid can drastically alter single-particle charge and spin 
transport. We show here that collective optical excitations arising from Coulomb interactions also 
exhibit unique signatures of Berry curvature and quantum geometric tensor. A non-zero Berry 
curvature mixes and lifts the degeneracy of Z 7 ^ 0 states, leading to a time-reversal-symmetric 
analog of the orbital Zeeman effect. The quantum geometric tensor, on the other hand, leads 
to Z-dependent shifts of exciton states that is analogous to the Lamb shift. Our results provide 
an explanation of the non-hydrogenic exciton spectrum recently calculated for transition metal 
dichalcogenides. Numerically, we find a Berry curvature induced splitting of ~ 10 meV between the 
2 pa; ± i2py states of WSe 2 . 


Introduction — An exciton, comprising of a bound 
electron-hole pair, is an elementary optical excitation of 
a semiconductor. In most semiconductors, large dielec¬ 
tric constant and small effective masses of charge carri¬ 
ers result in a weakly bound exciton with a Bohr-radius 
much larger than the lattice constant. Such excitons 
are termed as Wannier-Mott excitons and play a central 
role in understanding the optical response of a number 
of condensed-matter systems mm- Although the charge 
carriers in a semiconductor are described by Bloch-waves, 
excitonic spectra of most semiconductors resemble the 
hydrogenic series consisting of Rydberg series [3]. At a 
first glance, this is surprising since the wave functions 
of electrons and holes in a solid obey crystal symmetry. 
For example, discrete translational symmetry in a crystal 
implies that single-particle band dispersions of electrons 
and holes are periodic in crystal momentum k, necessar¬ 
ily making their kinetic energy non-parabolic. Neverthe¬ 
less, in a semiconductor with a direct band gap at the 
F-point (k = 0) and large Bohr radius such that exciton 
is made from k-states close to band-minima, full symme¬ 
try of vacuum is restored and a hydrogenic spectrum for 
excitons is obtained. 

It was recently reported that in monolayers of semicon¬ 
ducting transition metal dichalcogenides (TMDs) such as 
M 0 S 2 or WSe 2 , excitonic spectra deviates strongly from 
that of a hydrogen atom. In particular, states with iden¬ 
tical principal quantum number, such as 2s and 2p are 
not degenerate [US]- Moreover, the degeneracy of 2p 
states was also found to be lifted in tight-binding calcu¬ 
lations m- These exciting predictions were attributed 
to non-local dielectric screening in TMDs where exciton 
Bohr radius of the Is exciton is comparable to the mono- 
layer thickness m- In fact, as a consequence of the 
tightly bound nature of Wannier-Mott excitons in this 
material system, the exciton comprises of electron-hole 
states which are spread over a large momentum range 
where a parabolic description of the dispersion is not 
valid. Since TMD excitons are composed of electron- 
hole states around ±iF-points where Berry curvature is 


finite, it is natural to ask if Bloch-band geometry can 
alter the excitonic spectrum. While the role of Berry 
curvature in determining transport properties of non¬ 
interacting Bloch electrons is well established [HIS] , its 
role in modifying Coulomb attraction leading to bound- 
state formation has not been explicitly analyzed. 

In this Letter, we analyze the excitonic signatures of 
the two geometric invariants of Bloch bands - Berry cur¬ 
vature and quantum geometric tensor (QGT). We show 
that the Berry curvature is primarily responsible for a 
reciprocal-space orbital Zeeman effect leading to a finite 
splitting of the 2py;±i2py states. QGT on the other hand, 
contributes to the 2s-2p splitting, that is similar to the 
Lamb shift. Our findings apply to Wannier-Mott excitons 
in general, and are particularly relevant for TMD exci¬ 
tons where optical spectroscopy can directly probe the 
above-mentioned signatures of the Bloch-band geometry. 

Exciton problem in momentum space — The exci¬ 
ton motion can be decomposed into the relative mo¬ 
tion of electron and hole giving rise to hydrogen-like 
bound states, and the center-of-mass momentum (Kcm) 
resulting in an excitonic dispersion periodic in recip¬ 
rocal lattice vector. As light emission and absorp¬ 
tion takes places around Kqm = 0 due to negligi¬ 
ble momentum of photons, in the following we will 
only consider excitons with zero center-of-mass momen¬ 
tum. Without loss of generality, we restrict ourselves 
to two-dimensional excitons for the rest of the discus¬ 
sion. The exciton wave function can be expressed as 
- ■^ith Cc (c„) be¬ 

ing the annihilation operator for an electron in the con¬ 
duction (valence) band, and |0) as the semiconductor 
vacuum with no excitations. Amplitudes Ai,(k)s sat¬ 
isfy the eigenvalue equation + R|k')A,y(k') = 

AA,,(k) where U = Ek^k^lkC^.k + Ek'^’kCc.kCc.k is 
the single-particle Hamiltonian of the two-band semicon¬ 
ductor, and V = f (i^rcZ^r'dtl'(r)4'l(r') jjjri 77 |'I'(i’)di'(r') is 
the Coulomb interaction. Field operators are defined as 
" = ^E„=«,cEk«n,k(r)e**" '’c„,k, where M„.k(i-) 
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are the Bloch functions of the band with index n = c,v 
and S is the quantization area. 

The matrix elements are given by 


= |k\'| (^c.k|Mc,k')('»«,k'|'»«,k), which is the 

direct part of Coulomb interaction. The long-range 
part of exchange interaction of V vanishes for excitons 
with Kcm = 0. We neglect the short-range part of the 
exchange interaction since its magnitude is much smaller 
than the direct terms we analyze. 

The eigenvalue equation for the exciton then becomes, 

(Ak - E^)A^{k) - 1 ^ ^,4, ^,(k') = 0, 

( 1 ) 

where the index v = Is, 2s, 2p, etc in analogy to hy- 
drogenic orbitals, Ak = and Bloch overlaps 

®k k' = {un,'k\un.'k')■ When the Bloch overlaps are unity, 
and the dispersion Ak is parabolic, we recover the two- 
dimensional (2D) hydrogen atom solution. 

Bloch part: General arguments — It is clear from Eq. 

g that Bloch overlaps sjj enforce the symmetry of 
; crystal on Coulomb attraction between electron-hole 
pair. The radial symmetry of hydrogen atom is lowered 
when sj^ j., deviate significantly from unity thereby nec¬ 
essarily making the excitonic spectrum non-hydrogenic. 
This should be contrasted with a deviation arising from 
the non-parabolicity of single-particle states. Note that 
Bloch overlaps as defined above are gauge-dependent un¬ 
like the eigen-energies They need not change with 
band dispersion but depend on how the Bloch functions 
are arranged in k-space or the geometry of Bloch bands 
Un,\^ over the Brillouin zone. Guided by this observation, 
below we express sJJ j., in terms of geometric invariants of 
Bloch bands such as Berry curvature. We first note that 
the Bloch overlap for k' ^ k -f dk can be expressed 
as - 

Sk.k+dk = 1 +{Un{k)\dk,\Un{k))dki (2) 

-I- ^{Un{k)\dkidkj\un{k))dkidkj + ... 

The first order term in dk is the Berry connection 
iAi which is related to Berry curvature as r2(k) = 
V X ^(k) m Eq. ([^ is not gauge-invariant in that 
it does not transform as a tensor under the transfor¬ 
mation u„(k) = However, if one chooses a 

closed path in k-space, arbitrary phases which arise un¬ 
der the 17(1) gauge-transformation, mutually cancel each 
other to give gauge-invariant quantities. Indeed, only 
such closed-loop terms appear in the characteristic poly¬ 
nomial of the eigenvalue problem in Eq. Q. It can be 
shown that up to second order in dk (see Appendix A.) - 

(uk, |wk.)... (Mk„_, |uk„) (3) 


:|H|k') = d, 

^ / d2rd^/M*^(r)Mc,k^ 


k.k' 


and 
') 


£|r-r'| 


(k|E|k') = 

<,k'(i’')w«.k(r') 


for a closed path such that Wkjy = Wki- The first expo¬ 
nential on RHS is nothing but the Berry phase of the 
closed path while the second term in the exponential is 
the squared “length” of the path defined in terms of the 
quantum geometric tensor = Re [{dk^ u(k)|afc^.M(k))]- 
AiAj. QGT, also referred to as Fubini-Study metric, is 
a gauge-invariant quantity corresponding to the second 
order derivative in Eq. ([^ which measures infinitesimal 
distance between Bloch states parametrized by k m- 
Thus, the gauge-invariant Bloch overlaps over a closed 
loop can be expressed solely in terms of geometric quan¬ 
tities characterizing the Bloch bands. It is instructive 
to study the effect of Bloch overlaps perturbatively for 
which we consider an infinitesimal loop, when all ks are 
close to each other. Eq. (§ then becomes - 

1 -I- in • dSfc - ^g^jdkidkj -|- ..., (4) 

where we have expressed ^ A • dk in terms of Berry cur¬ 
vature as f fl ■ dSk using Stoke’s theorem. It should be 
noted that the imaginary part of the above expression 
is proportional to Berry curvature which is an antisym¬ 
metric quantity whereas the real part is proportional to 
QGT, a symmetric quantity. When calculating Bloch 
overlaps for conduction and valence bands, the difference 
(sum) of the Berry curvatures (QGTs) of the two bands 
appears. 

Illustrative model: Gapped graphene — Having es¬ 
tablished a connection between Bloch overlaps and ge¬ 
ometry of Bloch bands, we now illustrate how Berry cur¬ 
vature and QGT effect the exciton spectrum using a toy- 
model. We choose a simple two-band model of graphene 
with a band gap at iAT-points for which Bloch over¬ 
laps sJJ j./ can be analytically obtained. The Hamilto¬ 
nian can be written in Pauli basis of the two bands as 
Hik) = {atkx, atTyky, Ag) where a is the lattice constant, 
t denotes the hopping energy, and = ±I is the val¬ 
ley index. The presence of inversion symmetry breaking 
band gap Ag results in equal but opposite Berry curva¬ 
ture near the valleys at idf-points [nmi]. Furthermore, 
the magnitude of Ag is chosen to be large enough such 
that the magnitude of Berry curvature at the iAT-points 
|Dg| = a^t^/Ag is small and can be taken to be constant 
in the k region where the exciton wave function extends 
i.e., Dg|kp ^ 1 for |k| S dfc ^ l/as around the iAT- 
points. 

Under these assumptions, 

■S^k'SkCk = 1 + ^*T-„(k' X k) — -|k — k'l^^ -I- ..., 

(5) 

up to first order in |Hg|. Thus, |Dg| serves as a small 
parameter for perturbative treatment of the Bloch part. 
Note that Eq. (§ is analogous to Eq. Q for the present 
model. In particular, the imaginary, antisymmetric part 
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(opposite in the two valleys) corresponds to the Berry 
curvature while the real symmetric part to QGT. In the 
present model, both Berry curvature and QGT are pro¬ 
portional to Iflol = Plugging Eq. (|^ in Q, one 

can write the exciton Hamiltonian of Eq. Q perturba- 
tively as i7ex = , where is the Hamiltonian 

without the Bloch part describing the 2D hydrogen atom 
while is the perturbative term due to Bloch overlaps 
which reads as - 




1^01 27re^ / 

2 Ss \2 


k'l 


- ITv 


(k' X k)\ 


( 6 ) 


We decompose into a real, symmetric part, Vg and 
an imaginary, antisymmetric part, V^g, making Her- 
mitian, as expected. 

In order to determine the new eigenvalues and eigen¬ 
states of the perturbed Hamiltonian, one needs calculate 
the matrix elements of Vj in the basis of 2D hydrogenic 
levels. We first consider the effect of on 2px and 2py 
states which are degenerate in 2D hydrogen atom. Due to 
the symmetry of V^g under the exchange of k and k', its 
diagonal matrix elements vanish, while the off-diagonal 
matrix elements are finite. Likewise, Vg has only diago¬ 
nal matrix elements which are non-zero. 

The new eigenstates are obtained by diagonalizing the 
following matrix - 

[Eg + \Vi (V'2p. \Vig\^2pJ \ 

(7) 

V {i’2pJVig\i;2pJ Eg + |E/ \^2pgJ 


It is clear that the off-diagonal matrix elements arising 
from the antisymmetric part of the will mix the 2px 
and 2py states into symmetric and antisymmetric combi¬ 
nations \ip 2 p±) = {\^ 2 pJ ± i\i’ 2 py)) and causing them 

to split. On the other hand, the diagonal matrix elements 
due to the symmetric part of will cause rigid shift in 
energies of \'ip 2 p±) by the same amount. 

The above scenario is an analog of orbital Zeeman ef¬ 
fect where degeneracy of 2p states of a hydrogen atom 
is lifted in the presence of a constant magnetic field due 
to coupling to angular momentum I, an antisymmetric 
quantity. Berry curvature is the momentum space ana¬ 
log of magnetic field and the splitting here can be thought 
of as a “momentum-space orbital Zeeman effect”. In gen¬ 
eral, the degeneracy of alH 7 ^ 0 states of 2D hydrogen 
atom will be lifted due to Berry curvature as in the case 
of orbital Zeeman effect. Due to the appearance of val¬ 
ley index in E 45 , the splitting is opposite in the two 
valleys as required by time-reversal symmetry. Erom Eq. 
Q, we can conclude that when Berry curvature is small 
but not constant, it is the Berry flux through the exciton 
wave function which determines the splitting. We note 
the similarity of our findings with that of Ref. m where 
excitons on the surface of topological insulators with ex¬ 
plicitly broken time-reversal symmetry were considered. 


Under the assumption of small, constant flo and k close 
to RT-point, we numerically evaluate the matrix elements 
in Eq. 0 to be {'ip 2 pJVig\'ijj 2 py) = -it„|Do|c/2 and 
{'4’2p:^\Vs\^2p:^) ^ ~l^o|c/4, where c > 0 is a constant. 
The energies of 2p± states are obtained to be E 2 p^ = 
E 2 p — |no|c/4 ± |Ho|c/2. The splitting between the 2p 
states is then A 2 p oc r„|r2o|. Eor s-states, only the sym¬ 
metric part Vg survives due to symmetry reasons, lead¬ 
ing to a blue-shift in energy. We find that the shift of 2s 
states is ~ |r2o|c/4, making it almost degenerate with 2 p+ 
state under our assumptions. Figurej^ shows a schematic 
energy-level diagram of the model with and without the 
Bloch perturbation. Upon setting the Bloch overlap to 
be unity, we find that the 2p states remain degenerate 
resembling unmixed 2py; and 2py orbitals thereby con¬ 
firming the role of Bloch overlap in splitting. 

One can also consider a two-band model with iden¬ 
tically zero Berry curvature to confirm that the mixing 
of 2px and 2py states stems from Berry curvature. If 
{hx (kx, ky ), hy (kx, ky ), Ao) is the Hamiltonian in Pauli 
matrix basis, then choosing hx{kx,ky) oc hy{kx^ky) 
gives vanishing Berry curvature as it is proportional to 
dhx X dhy. Indeed, in numerical calculations we find that 
the 2p wave functions for such a model remain unmixed. 

Next, we comment on the role of QGT in determining 
the exciton spectra. QGT can be rewritten as - 

Qtj = R-e [{dkMdkju)] - {dkME){u\dk^u) (8) 
= {X^Xy) - {X,){Xy), 

where, the operator X is the generator of translation in 
the k-space. QGT then measures the quantum fluctua¬ 
tions of X M- In the present case, X corresponds to a 
spread in “relative position” of electron and hole. This is 
not unlike the case of Lamb-shift in hydrogen atom where 
vacuum fluctuations of the electromagnetic field smear 
electron’s position, in turn changing its potential energy 
and resulting in a measurable shift of its energy [TB]. In 
other words, QGT can be thought of as “vacuum fluc¬ 
tuations” of the effective U{1) gauge field present in the 
k-space, and the corresponding shifts as an analog of the 
Lamb shift. Even when Berry curvature vanishes iden¬ 
tically, the effect of QGT must still remain. It is note¬ 
worthy that the magnitude of this Lamb-like shift can be 
relatively large compared to the hydrogenic Lamb shift; 
this observation can be thought of as a consequence of 
the fact that the effective fine structure constant of the 
gapped graphene model, a = /hv is on the the order of 
unity. Thus, we identify another physically relevant con¬ 
sequence of QGT which has been previously shown to be 
behind such varied phenomena as electric polarization 
in insulators |19j . current noise |13] . and even quantum 
phase transitions nnun] 

TMD excitons — In the following, we investigate how 
the predictions of the preceding discussion apply to TMD 
excitons. We first note that the experimentally observed 
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FIG. 1. (Color online) A sketch of energy levels of the 
gapped-graphene model under assumptions of small, constant 
Berry curvature (flo) and large exciton Bohr-radius. In ab¬ 
sence of Bloch term in the numerical calculations, n = 2 states 
remain degenerate like hydrogen atom. Upon the inclusion of 
Bloch term, 2p states mix and split due to flo much like the 
orbital Zeeman effect. The shift in energy levels arises from 
quantum geometric tensor {gij) and is analogous to Lamb 
shift in hydrogen atom. 

deviation from a hydrogenic series arises due to a com¬ 
bination of non-local dielectric screening, non-parabolic 
dispersion of bands, and the effect of the Bloch part. 
The assumption of constant and small Berry curvature 
no longer holds due to the large spread of exciton wave 
function in k-space. In addition, unlike other material 
system such as GaAs where the lowest energy exciton is 
made from electron-hole states near the T-point (k= 0 ), 
lowest energy excitons in TMDs are made from ±Ff-point 
electron-hole pairs where there is non-zero Berry curva¬ 
ture. Thus, we expect a mixing and splitting of 2p-states 
in addition to energy shifts due to QGT. 

We solve for the eigenvalue problem in Eq. ([^ using a 
three-band, next-nearest neighbor model, for M 0 S 2 and 
WSe 2 , which captures the band dispersion of conduction 
and valence band throughout the Brillouin zone [55] . We 
discretize the Brillouin zone into a grid of 136 x 136 k- 
points and assume an air-suspended sample with a non¬ 
local dielectric screening length r* ^ 15 A, corresponding 
to a binding energy of ~ 400 meV for the Is-state which, 
is estimated to be in the range of 300 - 700 meV [51171155]. 
The Coulomb part of the Hamiltonian is calculated such 
that |k—k'l is always restricted to the first Brillouin zone. 

Figure]^ (a)-(b) show the 2p wave functions for M 0 S 2 
exciton with and without the Bloch part. Just like in the 
case of gapped graphene, 2p states remain unmixed with¬ 
out the inclusion of Bloch part in calculations. With the 
Bloch term included, we obtain a 2p splitting ^ 10 meV 
14 meV) for M 0 S 2 (WSe 2 ) which is consistent with 
the recently reported splitting of ^ 25 meV in M 0 S 2 m- 
Given the relatively large splitting, it should be possi¬ 
ble to detect it experimentally using two-photon, polar¬ 
ization resolved, optical spectroscopy involving a near- 
infrared laser slightly detuned from the Is exciton and 
a mid-infrared laser as shown in Fig. Sb). When the 



FIG. 2. (Golor online) Galculated squared amplitude of 2p 
wave functions in reciprocal space for M 0 S 2 in (a) absence 
and (b) presence of Bloch overlap. The Brillouin zone is cho¬ 
sen such that F-points lie on the four vertices. The exciton 
wave function extends near ±K points. Without the Bloch 
perturbation, the wave functions are degenerate and resemble 
2px and 2py wave functions. Bloch perturbation mixes the 
states which resemble 2p± states and have an energy split¬ 
ting of about ~ 10 meV. (c) A sketch of energy levels with 
a scheme to optically determine the predicted splitting of 2p- 
states using two-photon resonance spectroscopy. 

two lasers are co-circularly (counter-circularly) polarized, 
lower (higher) energy 2p state will be excited in a two- 
photon resonance leading to an enhancement in photo¬ 
luminescence from the Is state. We note that even in 
the absence of Bloch part, trigonal warping of the disper¬ 
sion can lead to a degeneracy lifting between 2px and 2py 
states, however, they remain unmixed and the magnitude 
of splitting is much smaller than the Berry curvature in¬ 
duced splitting. We also note that the level diagram in 
Fig. I^b) is strictly true for Kqm = 0, as the long-range 
part of exchange interaction will mix the states in oppo¬ 
site valleys for hnite Kqm [55] . 

Conclusions — The central role played by the Berry 
curvature in determining the transport properties of non¬ 
interacting Bloch electrons, leading to anomalous, val¬ 
ley and spin Hall effects, is well established. Our re¬ 
sults on the other hand, unequivocally demonstrate that 
the spectrum of Coulomb-correlated two-particle bound 
states exhibit observable signatures of Berry curvature 
and QGT. This motivates the question whether many- 
body optical excitations of semiconductors, such as trions 
in a 2D electron system, are also influenced by effective 
gauge-helds in solids arising due to non-trivial geometry 
of Bloch-bands. 
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APPENDIX A. 

Consider a (A^-l)-point discretization of a 
smooth curve C in k-space parametrized by points 


ki, k 2 ,..., kjv-i. We want to show that Eq. (3) in main 
text holds up to second order in dk. As we wish to 
express the Bloch overlap, sc as an integral over C, we 
take the logarithm of LHS to get, 


In(sc) = y^ln(ukjMk.+i)- 


(9) 


In the continuum limit. 


ln(uk|uk+dk) = \n{l+{u\di\u)dki+-{u\didj\u)dkidkj +...), 

where partial derivatives are taken in k-space and sum¬ 
mation over indices is assumed. Using ln(l J- cc) = 

X — . with X = {u\di\u)dki + ^{u\didj\u)dkidkj, 

we get up to second order, 

ln(uk|uk+dk) = {u\di\u)dki + ^{u\didj\u)dkidkj (11) 
-l^{u\di\u) {u\dj\u)dkidkj + ...), 


With the above identity, Sc becomes, 


Sc --exp 

® {u\di\u)dki 


Jc 


( 12 ) 


X exp 


- j) {{u\didj\u) — {u\di\u){u\dj\u))dkidkj ... 


Using Ai = —i{u\di\u) and Aidki = f fldki A dkj we 
get. 


sc exp 


Qdfci A dkj exp - j) {{u\didj\u) + AiAj) dkidkj 

(13) 

Under local U{1) transformation Uk = 

{u\didj\u) + AiAj —>■ {u\didj\u) + AiAj -I- ididja. (14) 

As Sc is gauge-invariant under the transformation, when 
integral over closed C is preformed, the gauge-dependent 
term ididja must vanish. Evidently, the real part of the 
above equation is gauge-invariant as A^s are real. There¬ 
fore, the following should hold. 


sc =exp 

X exp 


n ■ dSc 


(15) 


^ j) (Re [{u\didj\u)] + AiAj) dkidkj 


Using Re [{u\didj\u)] = —Re [{diu\dju)\, we arrive at Eq. 
(3) of main text. 

















